Electric-field control of a hydrogenic donor's spin in a semiconductor 
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An AC electric field applied to a donor-bound electron in a semiconductor modulates the or- 
bital character of its wave function, which affects the electron's spin dynamics via the spin-orbit 
interaction. Numerical calculations of the spin dynamics of a hydrogenic donor (Si) embedded in 
GaAs, using a real-space multi-band k • p formalism, show the high symmetry of the hydrogenic 
donor state results in strongly nonlinear dependences of the electronic g tensor on applied fields. A 
nontrivial consequence is that the most rapid Rabi oscillations occur for electric fields modulated 
at a subharmonic of the Larmor frequency. 



Electronic ground states characterized by non-zero 
spin are attractive candidates for encoding quantum in- 
formation in a solid state system, and the use of elec- 
tric fields is an attractive method to address individual 
spins U] E]. When the ground state has a nonzero inte- 
ger spin it is possible to perform all needed spin oper- 
ations using electric fields alone [3l |J, whereas for ma- 
nipulation of spin- 1/2 electronic ground states at least 
a static applied magnetic field is required. Advances in 
focused ion beam single- ion implant at ion [5], as well as 
the use of a scanning tunneling microscope to implant 
a single ion with atom-scale precision [7], suggest that 
spin clusters and spin-based circuits consisting of large 
numbers of precisely positioned spins could be designed 
with near-atomic resolution. Proposals to control indi- 
vidual spin- 1/2 states in such an environment with lo- 
cal electric fields include changing the magnitude of the 
Lande g tensor to bring spins into resonance with an ex- 
tended AC magnetic field O |9l [TOl [11], moving spins in 
a fringe-field [12 or a hyperfine[TT gradient, modulating 
zero- field spin splittings [14, J^, and ^-tensor modulation 
resonance (^-TMR) [16] . ^-TMR uses the electric-field de- 
pendence of the Lande g tensor anisotropy to manipulate 
the spin, and so does not require microwave magnetic 
fields or nanoscale magnetic materials or nuclear polar- 
ization gradients. Although ^-TMR works by changing 
the orbital character of the wave function with an electric 
field, and thereby indirectly influencing the spin through 
the spin-orbit interaction, it does not require zero-field 
spin splittings (so ^-TMR could be performed in a sil- 
icon or diamond host). Predictions for quantum dots 
indicate control of the g tensor anisotropy can produce 
rapid Rabi oscillations and full Bloch-sphere control with 
a single vertical electric field [17]. 

The promising approach of ^-TMR has yet to be ex- 
plored for electrons bound to dopants. Shallow donors 
might seem a poor candidate for modulation of g ten- 
sor anisotropy, as they have cubically-symmetric g ten- 
sors in the absence of an electric field. Quantum dots, 
by contrast, have highly asymmetric g tensors that are 
very sensitive to applied electric fields^. However, we 
find the g tensors of electronic spins bound to donors de- 
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FIG. 1 (a) Proposed geometry for ^-tensor modulation, 
(b) The impurity potential is a Coulomb potential plus a 
central cell correction. The effect of applying an electric field 
on the zone center band energies is shown. 



pend nonlinear ly on applied electric and magnetic fields, 
and thus substantial g tensor anisotropy and rapid spin 
manipulation can be achieved for a hydrogenic donor 
state. As the dominant electric-field dependence is non- 
linear, the most rapid Rabi oscillations are found at un- 
expected frequencies — subharmonics of the Larmor fre- 
quency rather than the fundamental — permitting rapid 
spin manipulation using AC electric fields with frequen- 
cies far below the Larmor frequency. Furthermore, the g 
tensors of quantum dots are very sensitive to dot shape 
and composition [18 and thus each quantum dot will have 
different resonance frequencies for ^-TMR. Donor wave 
functions and g tensors will, however, each be reliably 
the same. 

These hydrogenic states have other attractive features 
for spin clusters or spin devices; they possess the biggest 
radii of any ionic bound states in the solid, with Bohr 
radii of the order of 10 nm in GaAs. Thus the spin-spin 
coupling between states would be easier to control than 
for deep levels whose interaction strength changes sub- 
stantially on the atomic scale [7 . Our treatment focuses 
on the substitutional silicon donor in gallium arsenide, 
Sica, as it is one of the best understood semiconduc- 
tor point defects and is well described by the hydrogenic 
model. We expect that similar results are possible for 
a shallow donor in silicon, although the details may be 
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complicated by the presence of multiple valleys in the 
conduction band. 

The geometry of ^-TMR for a single electron spin 
bound to a donor is shown in Fig. [ija). A static mag- 
netic field along with a gated time varying electric field 
is applied to the crystal containing the Sica donor. We 
considered all orientations of the field and found that the 
most rapid Rabi oscillations occur when the magnetic 
field is applied at an angle = 45^ to the electric field, 
which is the configuration shown in Fig. [TJa). 

Although many properties of shallow impurities (such 
as the energy spectrum) can be treated to an excellent 
approximation by two-band effective mass theory [19 , g- 
tensor calculations require a multi-band treatment as 
the coupling among multiple bands needs to be con- 
sidered, and the spin-orbit interaction must be treated 
accurately!^. Moreover, the electric field breaks the 
spherical symmetry of the impurity site. These complex- 
ities are best handled numerically. 

Our calculations of ^-TMR for Sica donors were car- 
ried out using 8-band k- p theory [21 in the envelope ap- 
proximation using finite differences on a real space grid 
[IHl [22] [23l |24]. Material parameters were taken from 
Ref . [25l assuming T = 0. The potential of the hydrogenic 
impurity. 



K(r) = 



47rer 



C5(r-ro), 



(1) 



is the sum of a screened Coulomb potential and a 
delta function potential corresponding to the central-cell- 
correction (CCC). The CCC arises due to the differing 
chemical nature of various impurities. For our calcu- 
lations V^(r) is non-zero only on a single grid site, as 
shown in Fig. [ijb) along with the potential due to the 
applied electric field. The CCC is found by adjusting C 
until the calculated binding energy for the Is donor state 
matches experiment. The Lande g tensor for the impu- 
rity ground state was then obtained from the calculated 
Zeeman splitting of the Is level in a uniform magnetic 
field. 

Fig.|2]shows ^[001] 7 the tensor component for the Is im- 
purity state as a function of collinear magnetic and elec- 
tric fields. Increasing the electric field increases the rela- 
tive change in ^, whereas increasing the magnetic field de- 
creases the relative change in g. Appreciable changes in g 
are seen even at modest magnetic fields, which is encour- 
aging for manipulating the donor atom's spin. The im- 
purity 5'[ooi]'s depend nonlinearly on the magnetic field, 
as shown in Fig. [2|b) . This behavior is unlike that seen 
in small QDs such as treated in Ref. JX foi" which the 
g tensor is nearly independent of the applied magnetic 
field. 

The competing effects of B and E on ^[ooi] can be 
understood by examining the donor electron's wavefunc- 
tion, shown in Fig. [Sj As the magnetic field is increased 
in the [001] direction [from Fig. |3ja) to (b)], the cy- 
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FIG. 2 Normalized donor ^'[ooi] values as a function of (a) 
E[ooi] and (b) B[ooi]- Insets show unnormalized donor ^'[ooi] 
values. The full range of the x-axis of the insets is the same 
as that of their respective (normalized) plots. 




FIG. 3 Calculated wavefunctions for (a) 

^[001] = ^[001] = 0, (b) ^[001]= 0, B[ooi]= 4 T (c) 

^[001] = 0-15 mV/nm, B[ooi]= 0- Contours outline selected 

amplitudes as a guide to the eye. 



clotron radius decreases, contracting the extent of the 
wave function in the direction transverse to B. Simi- 
larly, the opposite effect is evident when the electric field 
is increased [Fig. [sj^c)], which allows the impurity wave- 
function to spread into a region with lower overall po- 
tential. This decreases the confinement for the donor 
electron and thereby increases |^|[18 . This effect is more 
prominent for a smaller magnetic field. 

The g tensor components were calculated for various 
directions of B with E applied along [001], as shown in 
Fig. [4j Note that dg/dE decreases with increasing B. 
The variation in dg/dE as a function of B is greater when 
E ± B. However at an intermediate B{^ 2T), dg/dE 
is identical in all directions. These results imply that 
an electric field induces a g tensor anisotropy oriented 
relative to E, which makes it possible to modulate the 
g tensor using an alternating electric field in addition to 
the static electric and magnetic fields. 

We next solve for the donor atom's spin dynamics 
by explicitly integrating the time-dependent Schrodinger 
equation. The nonlinear nature of g complicates a quanti- 
tative treatment within the rotating- wave approximation. 
The directional dependence of g (Fig. [4| can be used to 
obtain an analytical form of the g tensor by fitting each 



tensor component to the expression y ^ an{^)E 
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time-dependent g tensor can then be constructed for the 
time- varying electric field E(t) = Edc + Eac sin(a;t). The 
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FIG. 4 ^ as a function of ^[ooi] and B applied in various 
directions (a) [100] (b) [101] (c) [001]. g is also calculated 
for various magnetic field strengths. 

maximum amplitude of E{t) is always held constant at 
0.2 mV/nm, so as not to exceed the breakdown field of 
the GaAs host. The spin dynamics of the donor atom 
can then be calculated using the effective time-dependent 
Hamiltonian, 

H{t) = ^a-g{t)-B (2) 

where is the Bohr magneton. 

As the Hamiltonian is explicitly time dependent, the 
state of a spinor, Sj (where j =i,T) at time t can be 
obtained by evolving Sj (t = 0) forward in time in n steps 
of At = t/n <^ 1/200 as follows, 




where T is the time-ordering operator. For sufficiently 
small At this is equivalent to 

|5,(t)) = texp^*(^^i^)|5,(0)), (4) 

The time dependent probability of making a spin-flip 
transition is |(S'|(0)|S'|(t))p. Rabi oscillations are ob- 
tained when spin flip transitions are made resonantly 
(i.e. \{S^{0)\Si{t))\'^^^ = 1). Resonant spin flip transi- 
tions are usually made when E{t) is driven at the Larmor 
frequency Ql- However in case of the hydrogenic impu- 
rity system considered here, the donor electrons spin can 
be resonantly flipped at any sub-harmonic of the Larmor 
frequency: Ql/N^ where N is an integer. This is illus- 
trated in Fig. |5|a), where the peak spin-flip transition 
probabilities are shown as a function of the driving En- 
field frequency co. Multiple resonance lines are apparent, 
located at and its sub-harmonics. This unusual be- 
havior arises from the highly nonlinear dependence of g 
on the applied electric field (Fig. [4|. For sub-harmonics 
higher than N=2, the Rabi frequencies Qr are lower than 
those at A/" < 2 and hence are not considered further for 
spin manipulation. The largest Qr can be achieved by 
driving E at the second sub-harmonic {N = 2) of Ql- 



Due to the smaller DC component of the electric field 
the Rabi oscillations are less rapid at than at its 
second sub-harmonic. The resonance lines in Fig. [sja) 
at CO = VLl/N^ have a full width at half maximum of 
Auo = 2nR/N. 

The Rabi frequencies are calculated next as a function 
of Edc and and are shown in Fig. jsjb) with the electric 
field driven at l^^- For all ^, and uj = Qr is largest 
when the AC and DC components of the electric field are 
equal. If the electric field is driven at 1^l/2, however, as 
shown in Fig.jsjc), then Qr is largest if E^c = 0. In both 
Figs.jsjb) and (c), the optimal angle of the magnetic field 
to the electric field is ^ = 45^. Although the maximum 
Qr in Figs.|5|b) and (c) are identical, driving £^ at ^7^/2 
offers two key advantages. When the peak value of E is 
close to the breakdown of the host crystal, a pure AC 
field with an adjustable duty-cycle is much less likely 
to ionize the donor electron, as the carriers can recover 
during a thermal relaxation time. This allows for higher 
driving fields, which result in higher Qr. It also may be 
experimentally more feasible to resonantly flip the spin 
at the lower frequency of the subharmonic than the 

fundamental Ql- 

Fig. jsjd) shows Qr as a function of B for = 45^ 
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FIG. 5 Spin dynamics of the donor atom as a function of 
various parameters. Eac + Edc = 0.2 mV/nm and is [001] 
oriented. is the angle between B and E. (a) Peak spin- flip 
transition amplitudes as a function of E"s driving frequency, 
for Eac/ Edc = 9 and = 45°. Resonant transitions appear 
at sub-harmonics of the Larmor frequency Ql- (be) Rabi 
frequency Qr as a function of Edc and for: (b) E driven at 
Ql, is maximum at ^ = 45° and Edc — 0.1 mV/nm. (c) 
E driven at Ql/2, Qr is maximum at ^ = 45° and Edc = 0. 
(d) Qr as a function of B for optimal and Edc of (b) and 
(c). Note that above B — 2T, Q.r increases monotonically. 



4 



and E driven at Ql or Ql/"^- For magnetic fields greater 
than 2T, Qji increases monotonically, whereas below 2T 
Qr exhibits a non-monotonic behavior. This feature can 
be explained by Taylor-expanding the time-dependent 
Hamiltonian to first order in the rotating wave approxi- 
mation, 

Here the Larmor frequency is given by the time inde- 
pendent static precession vector, CIq = jub^ ■ 'B/h and 
the electron's spin dynamics in the rotating frame is de- 
scribed by the time-dependent spin precession vector, 
fti{t) = /j.BEac{dg/dE) • 'B/2h. Qi can be resolved 
into components that are parallel and perpendicular 
(r^x) to CIq. In the rotating frame, is equivalent to 
Qji (in the lab frame), as driving E at leads to spin 
precession about ft± or Rabi oscillations. As the ten- 
sor components dg/dE decrease with increasing B (see 
Fig. [4|, the magnitudes of B and dg/dE have opposing 
effects on l^i (and hence l^x). For B < YT the con- 
tribution from dg/dE dominates over B and hence the 
Rabi frequencies increase. For 1T< B < 2T the com- 
peting contributions of B and dg/dE make the g tensor 
increasingly isotropic and the Rabi frequencies smaller. 
At 5 2T the g tensor becomes isotropic and the Rabi 
frequency vanishes. Above B ^ the effects of a much 
larger magnetic field dominate and the spin flip times 
decrease monotonically. Two key inferences, consistent 
with other work on ^-TMR, can be drawn from this be- 
havior. For spintronic applications the highest magnetic 
field possible is desirable in order to generate the largest 
possible Rabi frequencies. Secondly, the amount of g- 
tensor anisotropy induced is crucial to achieving shorter 
spin-flip times, not the degree of change in ^ as a function 
ofE. 

We have proposed a scheme for achieving electric-field 
driven ^-tensor modulation resonance for a single shal- 
low donor impurity. Electric and magnetic field depen- 
dent g tensors were calculated for the Sica donor using 
8-band k • p theory on a real space grid. Varying E and 
B affects the confinement for the donor electron, which 
in turn alters its g tensor. In addition to the nonlinear 
E dependence, the g tensors are also highly nonlinear 
as a function of B. This is unlike the case for a QD, 
where g is essentially independent of B. A consequence 
of the nonlinear dependence of ^ on £^ is that spin-flip 
transitions can be made resonantly at any sub-harmonic 
of the Larmor frequency. Spin flip times were calculated 
exactly, using time evolution operators, and optimized 
for various parameters of interest. If E is driven at the 
second sub-harmonic of the Larmor frequency, then high 
frequency Rabi oscillations can be obtained without any 
DC component of E. This could be particularly use- 
ful in obtaining the largest Rabi frequencies for a given 
breakdown- field limit for the semiconductor host. 
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